Strangeness fluctuations and MEMO production at FAIR 



O 

o 
< 



I ' 

Oh- 



> 
O 



00 

O 



X 



Jan Steinheimer^ , Michael Mitrovski^'^, Tim Schuster^'^, Hannah Petersen^'^, Marcus Bleicher^ and Horst Stocker^'^'"^ 
^ Institut fiir Theoretische Physik, Johann Wolfgang Goethe- Universitat, 
Max-von-Laue-Str. 1, 60438 Frankfurt am Main, Germany 
^ Frankfurt Institute for Advanced Studies (FIAS), Johann Wolfgang Goethe- Universitat, 
Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany and 
^ GSI - Helmholtzzentrum fiir Schwerionenforschung mbH, Planckstr. 1, Darmstadt 

We apply a coupled transport-hydrodynamics model to discuss the production of multi-strange 
meta-stable objects in Pb+Pb reactions at the FAIR facility. In addition to making predictions 
for yields of these particles we are able to calculate particle dependent rapidity and momentum 
distributions. We argue that the FAIR energy regime is the optimal place to search for multi- 
strange baryonic object (due to the high baryon density, favoring a distillation of strangeness). 
Additionally, we show results for strangeness and baryon density fluctuations. Using the UrQMD 
model we calculate the strangeness separation in phase space which might lead to an enhanced 
production of MEMOs compared to models that assume global thermalization. 



Massive heavy-ion reactions provide an abundant 
source of strangeness. M ore t han 50 hyperons and about 
30 Anti-Kaons (i.e. -\-K^ carrying the strange quark) 
are produced in central collisions of lead nuclei at the 
CERN-SPS low energy program and before that at the 
AGS (see e.g. 01). In the near future, the Facility for 
Anti-proton and Ion Research (FAIR) will start to in- 
vestigate this energy regime closer with much higher lu- 
minosity and state-of-the-art detector technology. This 
opens the exciting perspective to explore the formation 
of composite objects with multiple units of strangeness 
so far unachievable with conventional methods. 

Exotic forms of deeply bound objects with strangeness 
have been proposed long ago (see ^!|) as collapsed states 
of matter, either consisting of baryons or quarks. For ex- 
ample the H di-baryon (a six quark state) was predicted 
by Jaffe 0. Later a multitude of bound di-baryon states 
with strangeness were proposed using quark potentials 
[1, Ol or the Skyrme model Q. However, the (non- 
)observation of multi-quark bags, e.g. strangelets and 
(strange) di-baryons is still one of the open problems of 
intermediate and high energy physics. Most noteworthy 
in this respect has been the hunt for the Pentaquark 
over the last 10 years, which re-stimulated this field and 
resulted in a reported observation at the CERN SPS 
accelerator 

The early theoretical models based on SU(3) and SU(6) 
symmetries [1, and on Regge theory [lol . [ll| suggest 
that di-baryons should exist. More recently, even QCD- 
inspired models predict di-baryons with strangeness S 
= 0, -1, and -2. The invariant masses range between 
2000 and 3000 MeV 0, [H, M, Q E H [13, [Hi • Un- 
fortunately, masses and widths of the expected 6-quark 
states differ considerably for these models. Nevertheless, 
most QCD-inspired models predict di-baryons and none 
seems to forbid them. 

On the conventional hadronic side, however, hypernu- 
clei are known to exist already for a long time [19, 23|. 
The double A hypernuclear events reported so far are 
closely related to the H di-baryon [2l| . Metastable exotic 



multi-hypernuclear objects (MEMOs) as well as purely 
hyperonic systems of A's and S's were introduced in 
[22, [2^ as the hadronic coun terp arts to multi-strange 
quark bags (strangelets) [13, HJ. Most recently, the 
Nijmegen soft-core potential was extended to the full 
baryon octet and bound states of EE, SS, and SS di- 
baryons were predicted [2^. For previous estimates of 
strangelet production and MEMO formation, the reader 
is referred to [13, El. 

A major uncertainty for the detection of such specu- 
lative states is their (meta) stability. Metastable exotic 
multi-hypernuclear objects (MEMOs), for example, con- 
sist of nucleons, A's, and S and are stabilised due to 
Pauli's principle, blocking the decay of the hyperons into 
nucleons. Only few investigations about the weak decay 
of di-baryons exist so far (see for a full discussion and 
new estimates for the weak nonleptonic decays of strange 
di-baryons) : In [2§| , the ff-di-baryon was found to decay 
dominantly hy H ^ T,~ -\- p for moderate binding ener- 
gies. While the (AA) bound state, which has exactly the 
same quantum numbers as the ff-di-baryon, was studied 
in 1131 ■ Here, the main non-mesonic channel was found to 
be (AA) ^ A-\-n. If the life time of the (AA) correlation 
or particle is not too long, the specific decay channels 
might be used to distinguish between both states. 

There are several searches in heavy-ion collisions for 
the i?-di-baryon j33,[3l,[33| and for long-lived strangelets 
[ssl. [3^ with high sensitivities, so far with no conclusive 
results. In pN collisions at the Fermilab however, the 
i/-di-baryon seems to be excluded over a wide range of 
Masses (2.194 < Mh < 2.231 GeV) and lifetimes (5 • 
10"^" to 1-10~^ sec) [31|. Hypernuclei have been detected 
most recently in heavy-ion reactions at the AGS by the 
E864 collaboration E^i. 



I. MEMO PRODUCTION RATES 

In this paper we study the production rate of multi- 
strange objects within the UrQMD model (v2.3) and a 
micro+macro hybrid approach to heavy ion reactions. 
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Similar to the RQMD model [H, [3§| which was em- 
ployed in UrQMD is a microscopic transport ap- 
proach based on the covariant propagation of constituent 
quarks and diquarks accompanied by mesonic and bary- 
onic degrees of freedom. It simulates multiple interac- 
tions of ingoing and newly produced particles, the exci- 
tation and fragmentation of color strings and the forma- 
tion and decay of hadronic resonances. At RHIC ener- 
gies, the treatment of sub-hadronic degrees of freedom 
is of major importance. In the UrQMD model, these 
degrees of freedom enter via the introduction of a for- 
mation time for hadrons produced in the fragmentation 
of strings [iol . l4ll . |43 | . The leading hadrons of the frag- 
menting strings contain the valence-quarks of the original 
excited hadron. In UrQMD they are allowed to interact 
even during their formation time, with a reduced cross 
section defined by the additive quark model, thus ac- 
countin g fo r the original valence quarks contained in that 
hadron p3.l4i|. 

For the microscopic+macroscopic calculation, the 
Ultra-relativistic Quantum Molecular Dynamics Model 
(UrQMD) is used to calculate the initial state of a heavy 
ion collision for a subsequent hydrodynamical evolution 
[i^ . H . [45| . This has been done to account for the 
non-equilibrium dynamics in the very early stage of the 
collision. In this configuration the effects of event-by- 
event fiuctuations of the initial state are naturally in- 
cluded. The coupling between the UrQMD initial state 
and the hydrodynamical evolution proceeds when the 
two Lorentz-contracted nuclei have passed through each 
other. 

_ 2R 

tstart — I (1) 

After the UrQMD initial stage, a full (3+1) dimen- 
sional ideal hydrodynamic evolution is performed using 
the SHASTA algorithm [IE, El. For the results pre- 
sented here an equation of state for a hadron-resonance 
gas without any phase transition is used [4^. The EoS 
includes all hadronic degrees of freedom with masses up 
to 2 GeV, which is consistent with the effective degrees 
of freedom present in the UrQMD model. One should 
note that we apply a purely hadronic EoS, for energy 
densities where a transition to the QGP is expected (see 
also [i^ for details on the model and comparison of 
extracted particle yields to data). Final particle (and 
MEMO) multiplicities are mainly sensitive on the de- 
grees of freedom at chemical freezeout which is refiected 
in the hadronic EoS. Dynamical observables such as 
momentum and rapidity spectra are more sensitive on 
the underlying dynamics. In addition, a phase transition 
could catalyse a strangeness destination process further 
enhancing MEMO production. However, studying the 
effects of a phase transition on MEMO production is left 
subject of future investigations. 

The hydrodynamic evolution is stopped, if the energy 
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TABLE I: Properties of all considered multibaryonic states 

density of all cells drops below five times the ground 
state energy density (i.e. ~ 730MeV/fm'^). This cri- 
terion corresponds to a T-/iB-configuration where the 
phase transition is expected - approximately T = 170 
MeV at fiB — 0. The hydrodynamic fields are mapped to 
particle degrees of freedom via the Cooper- Frye equation 
on an isochronous hyper-surface. 

j f{x,p)pfda^, with dCT^ = (dx3,0) (2) 

Here f{x,p) are the boosted Fermi or Bose distributions 
corresponding to the respective particle species. Inputs 
for these distributions are the masses and chemical po- 
tentials of the desired particles. For our calculation we 
assumed the mass of a MEMO to be the sum of the 
masses of all hadrons it is composed of. Similarly the 
total chemical potential is the sum of the constituents, 
and is composed of baryon and strange-quark chemical 
potentials and /is- 

The particle vector information is then transferred back 
to the UrQMD model, where rescatterings and the final 
decays are performed using the hadronic cascade. Using 
this parametrisation of the model one obtains a satisfac- 
tory description of data in a energy regime of 1 — 160A 
GeV. A more detailed description of the hybrid model in- 
cluding parameter tests and results for multiplicities and 
spectra can be found in fi^l- 

To calculate the multipHcities of MEMOS in the FAIR 
energy region, we employ the introduced hybrid approach 
to heavy ion collisions. Thus, the fiuctuating initial 
state produced in UrQMD, is coupled to a (3+1) di- 
mensional hydrodynamics evolution. When the energy 
density drops below 5eo(^ 730MeV/fm^) the freeze-out 
is performed and MEMOs and strangelets are produced 
according to the Cooper- Frye description |[2|). As dis- 
tinctive inputs for the distribution functions, the chem- 
ical potentials (/is,/^^) and masses of the MEMOs en- 
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FIG. 1: Multiplicities of various types of MEMOs and 
strangelets in central Pb+Pb reactions at -Eiab = 30A GeV 
from the hybrid approach. 



ter as discussed above. Final state interactions of these 
MEMOs are neglected for the present study. Table @ 
gives the properties of all multibaryonic states consid- 
ered in our analysis. They are the most promising and 
stable candidates. 

Fig. ^ provides the total multiplicities per degener- 
acy factor of various types of MEMOs and strangelets 
in central Pb+Pb reactions at i^iab — 3QA GeV. The 
yields obtained are in good comparison to the statistical 
model analysis [HOl, -which is describing strange cluster 
production at AGS energies. 

One should also note that -we assume particle production 
from a grand canonical ensemble for all beam energies. 
Because local, as well as global, thermal equilibration 
are assumptions not necessary justified in heavy ion 
collisions, a microcanonical description, combined with 
MEMO production by coalescence, has been proposed 
in Due to the restrictions of energy and momentum 
conservation, resulting in a phase space reduction for 
produced strange particles a (micro) canonical descrip- 
tion of the system strongly decreases strange particle 
yields (5l|,[53,[53. 

On the other hand, thermal models are able to repro- 
duce strange particle yields for beam energies above 
Eiab ~ GeV very well, and canonical corrections 
become negligible above these energies [H3|- 

Investigating strange-cluster production over a range 
of beam energies shows a distinct maximum in the yields 
of several multi strange objects. Fig. [2] displays the ex- 
citation function of the multipHcities of various MEMOs 
in central Pb+Pb reactions from the hybrid approach. 
The presented MEMO candidates are expected to pos- 



FIG. 2: (color online) Excitation functions of the multiplici- 
ties of various MEMOs in central Pb+Pb reactions from the 
hybrid approach. 

sess binding energies up to Eb/Ab ~ —22 MeV[25|. 
One easily observes that the upper FAIR energy region 
E\ah = 10 — 4QA GeV) is ideally placed for the search 
of exotic multi-strange baryon clusters. At lower ener- 
gies, the hyperon production cross section is too small, 
while at energies above FAIR, the expansion of the source 
and the small baryo-chemical potential suppress the for- 
mation of MEMOs and strangelets. 

Using the hybrid model enables us to also explore the 
phase space distribution of the produced particles. Fig. 
[3] shows the rapidity density of various MEMOs in cen- 
tral Pb+Pb reactions at E'lab — 30A GeV from the hy- 
brid approach. The production of baryon rich clusters 
is most pronounced in the high baryon density rapid- 
ity region. The rapidity distributions for MEMOs with 
a larger strangeness to baryon number fraction tend to 
look more gaussian like. 

Figure H] depicts the transverse momentum distribu- 
tion of various MEMOs at midrapidity in central Pb+Pb 
reactions at i?iab — 30^ GeV from the hybrid approach. 
The pt spectra are rather broad as compared to usual 
hadrons. This is due to the large boost the MEMOs 
acquire due to their large mass and the fact, that they 
are produced predominantly in the hottest regions of the 
expanding system. 



II. FLUCTUATIONS 

For the present study so far we have assumed global 
as well as local strangeness conservation. These assump- 
tions are common for models including thermal produc- 
tion of particles. In the following we explore if that as- 
sumption of local strangeness conservation is justified. 
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FIG. 3: (color online) Normalized rapidity density of various 
MEMOs in central Pb+Pb reactions at Siab = 30^1 GeV from 
the hybrid approach. 
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FIG. 4: (color online) Transverse momentum spectra at 
midrapidity {\y\ < 0.5) of various MEMOs in central Pb+Pb 
reactions at -Eiab = SOA GeV from the hybrid approach. 



especially in the FAIR energy regime. A relaxation of 
this assumption within the hybrid approach will require 
the expHcit propagation of the strangeness density (sim- 
ilar to the treatment of the baryon density). A second 
key ingredient will be the inclusion of an equation of state 
that can provide p{e, pB, Ps) with a finite ps- However, 
first we explore if such an extension might be necessary 
by applying the UrQMD model without an intermediate 
hydrodynamic phase. 

We start with an investigation of the strangeness 
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FIG. 5: (color online) Energy dependence of the strange 
quark over anti-strange quark (s/s) ratio for central 
Pb+Pb/ Au+Au reactions. Circles show the ratio at midra- 
pidity, while squares show the 47r values where the ratio is 
unity due to strangeness conservation. 



production and its distribution as a function of en- 
ergy. Fig. [5] depicts the energy dependence of the 
strange quark over anti-strange quark (s/s) ratio for cen- 
tral Pb+Pb/Au+Au reactions. The red circles present 
the strangeness to anti-strangeness ratio at midrapidity, 
while the blue squares show the 47r values where the ra- 
tio is unity due to strangeness conservation. One clearly 
observes that strangeness is not evenly distributed over 
rapidity, leading to an asymmetry between strange and 
anti-strange quarks on the level of 20% in the relevant 
energy regime. A similar kind of strangeness separation 
process has been predicted long ago within models cou- 
pling a hadron gas to a Quark-Gluon-Plasma state [13, 
|55| . Within these models the energy and particle num- 
ber balance in the mixed phase supports a 'distillation' 
process that enriches the QGP phase with strangeness 
and the hadronic phase with anti-strangeness. Within 
the present model, however, hadronic interactions are re- 
sponsible for the phase space separation of strangeness 
and anti-strangeness since no first order phase transition 
is present. Since both procedures separate strangeness 
in an equivalent way one can expect an even stronger 
strangeness separation if both effects are at work. 

Fig. [6] shows the rapidity dependence of the strange 
quark over anti-strange quark (s/s) ratio for central 
Pb+Pb/ Au+Au reactions at AGS {Eiab = 2 - llA GeV), 
SPS {Slab = 20 - 158^ GeV) and RHIC {^JsWn = 
19 - 200 GeV) energies. In the AGS and SPS energy 
regime, the (s/s) ratio is strongly rapidity dependent and 
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FIG. 6: (color online) Rapidity dependence of the strange 
quark over anti-strange quark (s/s) ratio for central 
Pb+Pb/Au+Au reactions at AGS {Eiab = 2 - llA GeV), 
SPS {Eiab = 20 - 158A GeV) and RHIC (^s^ = 19 - 200 
GeV) energies. 



FIG. 7: (color online) Transverse momentum dependence of 
the strange quark over anti-strange quark (s/s) ratio for cen- 
tral Pb+Pb/Au+Au reactions at AGS {Eiab = 2-11^ GeV), 
SPS {Eiab = 20 - 158yl GeV) and RHIC (^s^ = 19 - 200 
GeV) energies. 



has a pronounced peak above unity near midrapidity. At 
RHIC energies, the (s/s) ratio turns into a box shape 
as a function of rapidity with a plateau at unity indi- 
cating that strangeness is locally neutralized in rapidity. 
At the highest RHIC energies, the (s/s) ratio even turns 
slightly smaller than 1. Continuing this trend one vi^ould 
expect a clearly smaller than 1 (s/s) ratio at LHC ener- 
gies. In consequence, statistical model approaches (with 
the constraint of strangeness conservation at mid rapid- 
ity [i^l) are allowed to use midrapidity particle ratios as 
input for their calculations only at low RHIC energies. 
At lower, as well as higher energies, this procedure is not 
justified as strangeness neutralisation does not hold for 



the central rapidity region. Thermal calculations, using 
full phase space data as an input [H^ . [57| . and results 
from a thermal model including a core-corona scenario 
[H^ l. generally give better descriptions of strange parti- 
cle data, supporting the idea of dynamical strangeness 
separation, as proposed by this work. . 

Fig. [7] provides the transverse momentum dependence 
of the strange quark over anti-strange quark (s/s) ra- 
tio for central Pb+Pb/Au+Au reactions at AGS {Eiai, — 
2 - UA GeV), SPS {Eiab = 20 - 158A GeV) and RHIC 
{y/sNN = 19 — 200 GeV) energies. Here one observes 
a strong separation of strangeness in transverse direc- 
tion. With decreasing energy (increasing baryo-chemical 
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FIG. 8: Fluctuations of the strangeness fraction fs — Ps/pB 
in the central plane (x is in the impact parameter direction, y 
is transversal to the impact parameter and longitudinal direc- 
tion) for a single central Pb+Pb reaction at Eiab = SOA GeV. 
The colour coding indicates the local strangeness fraction, 
dark regions have more anti-strange than strange quarks. 



potential) the distribution of (anti-)strangeness becomes 
increasingly non-uniform in momentum space. The low 
momentum region is depleted of strange quarks, while the 
high pt region shows a strong enhancement of strange 
quarks compared to anti-strange quarks. This can be in- 
tuitively Hnked to the fact that (multi) strange baryons 
have a larger inverse slope than the Kaons for a given 
transverse velocity due to their larger masses.^ 

Next, we turn to the distribution and fluctuations of 
strangeness in coordinate space. Fig. [8] elucidates the 
fluctuations of the strangeness fraction fs — Ps/pB, with 
PB being the local baryon density and ps being the local 
net-strangeness density, in the central plane for a single 
central Pb+Pb reaction at Eiab = 30A GeV. Here x is in 
the impact parameter direction and y is transversal to the 
impact parameter and longitudinal direction. The distri- 
bution of the net-strangeness and baryon densities were 
obtained from the UrQMD model by means similar to 
creating the hydro initial state in the hybrid model. All 
hadrons and their baryon number and strangeness con- 
tent are represented by a Gaussian with a finite width of 
1 fm [isl . I49I . The plot is shown for the time when both 
nuclei have passed each other. The colour coding indi- 
cates the local strangeness fraction, white regions have 
more strange than anti-strange quarks, while dark and 
black regions show more anti-strange quarks. Locally 



strangeness lumps of 4 fm^ x Az appear both in positive 
and negative strangeness directions. As for the distribu- 
tion in momentum space discussed above, also the coordi- 
nate space distribution is largely non-uniform, although 
these spacial fluctuations occur only on an event-by-event 
basis. 

We have presented results for the thermal production 
of MEMOs in nucleus-nucleus collisions from a combined 
micro+macro approach. MultipHcities, rapidity and 
transverse momentum spectra are predicted for Pb+Pb 
interaction at Eiab = 5^1 GeV and Eiab = 30^ GeV. The 
presented excitation functions for various MEMO multi- 
plicities show a clear maximum at the upper FAIR en- 
ergy regime making this facility the ideal place to study 
the production of these exotic forms of multistrange ob- 
jects. Detector simulations have shown that the CBM 
experiment is well suited for the search of exotic multi- 
hypernuclear objects either by invariant mass reconstruc- 
tion of strange di-baryons or observing decay systematics 
(The very stable double negative {2S°,2S~} for exam- 
ple should have a characteristic decay in two negatively 
charged particles) [25]. 

Compared to many previous studies on MEMO and 
strangelet production, based on statistical models with 
global strangeness and baryon number conservation, the 
present approach indicates that the local strangeness 
density clumps strongly in coordinate space and that 
strangeness is unevenly distributed in momentum space. 
This mechanism does not require the production of a de- 
conflned state, and proflts from the non-equilibrium fea- 
tures present in the reaction. These fluctuations might 
lead to an enhancement of MEMO (and strangelet) pro- 
duction compared to previous calculations. 
Furthermore the net strangeness at midrapidity deviates 
from zero - not only on an event-by-event basis - indi- 
cating that the assumption of local strangeness neutral- 
isation is only justifled at the RHIC energy regime, but 
not at lower energies. Here it is therefore questionable 
if midrapidity particle ratios can be used as input for 
thermal particle multiplicity calculations. 
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^ If one assumes strangeness conservation as well as vanishing net 
baryon number at midrapidity (as is expected at very high en- 
ergies) then these distributions should be flat, as particles and 
their antiparticles are produced in equal numbers. If the (s/s) 
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